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Abstract 

Background: The detection of patfiogens at eariy stages of infection is a l<ey point for disease control in 
aquaculture. Tlierefore, accurate diagnostic procedures are a must. Real-time PCR has been a mainstay in 
diagnostics over the years due to its speed, specificity, sensitivity, reproducibility and throughput; as such, real-time 
PCR is a target for improvement. Nevertheless, to validate a novel diagnostic tool, correct setup of the assay, 
including proper endogenous controls to evaluate the quantity and quality of the samples and to detect possible 
sample degradation, is compulsory. This work aims to design a unique RT-qPCR assay for pathogen detection in the 
three salmonid species reared in Chile. The assay uses elongation factor 1 alpha as the single endogenous control, 
thus avoiding the need for multiple endogenous controls, as well as multiple validations and non-comparable 
quality control parameters. 

Results: The in vivo and in vitro analyses of samples from Salmo salar, Oncorliynclius myl<iss and Oncorhynclius 
l<isutcli showed that when primers were accurately selected to target conserved regions of the elongation factor 
1 alpha (ELFla) gene, a single novel RT-qPCR assay yielding similar and reproducible Ct values between the three 
species could be designed. The opposite occurred when an assay originally designed for Salmo salar was tested in 
samples from the two species of the genus Oncorhynchus. 

Conclusions: Here, we report the design and evaluation of an accurate trans-species RT-qPCR assay that uses the 
elongation factor 1 alpha (ELFla) gene as an endogenous control and is applicable for diagnostic purposes in 
samples obtained from the three salmonid species reared in Chile. 

Keywords: Elongation Factor 1 alpha, ELFla, Endogenous control, Fish disease. Fish virus, Real-time PCR, RT-qPCR 



Background transport, poor water quality and overstocking. These 
In nature, fish are exposed to a wide variety of microorgan- stressors, together with other weaknesses such as physio- 
isms of bacterial, viral or fungal origins, many of which are logical unbalance or nutritional deficiency, challenge 
capable of causing disease. When fish are reared under homeostasis and thus allow opportunistic infections to 
controlled conditions, the menace of disease-causing proceed [1]. As a result, diseases represent a major eco- 
agents increases significantly because confined fish are also nomic cost, especially when the impact is on fish of high 
exposed to a number of stressors, such as handling, commercial value [2] . Some of these pathogens are fastidi- 
ous and highly aggressive, causing high mortality rates, 
— — , , while others persist in fish, representing a potential danger 
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diseases that seriously threaten the sustainability of aqua- 
culture. One of the techniques that has consistently been 
transversal in specifically detecting sequences of target 
pathogens in humans, plants and animals, including fish, is 
real-time PGR because of its sensitivity, specificity, speed, 
throughput and reproducibility [3]. In the aquaculture en- 
vironment, this technique has been successfully used to de- 
tect, among other pathogens, viral hemorrhagic septicemia 
virus (VHSV) [4], infectious hematopoietic necrosis virus 
(IHNV) [5], salmon alphavirus (SAV) [6], Piscirickettsia 
salmonis [7-9], Renibacterium salmoninarum [10], infec- 
tious salmon anemia virus (ISAV) [11], and infectious pan- 
creatic necrosis virus (IPNV) [12]. 

Nevertheless, massive application of this technique re- 
quires correct setup of the assay, including the necessary 
controls, for unequivocal interpretation of the results. 
Among the most important of these controls is an en- 
dogenous reference control that is able to provide infor- 
mation not only about the amount of the detectable 
target nucleic acid in the sample to be analyzed but also 
about the sample integrity in order to accurately validate 
processes such as sampling, transport, and nucleic acid 
isolation, allowing degradation of the samples to be 
detected [11,13]. This is not a simple task to achieve, as 
determining the proper endogenous levels of a given 
marker molecule depends on a number of variables, in- 
cluding the target organ and the type of nucleic acid to be 
detected. In the case of targeting viral RNA, specifically 
for the diagnosis of the ISA virus in infected tissues of fish, 
the use of the elongation factor 1 alpha (ELFla) has been 
demonstrated to be a reliable reference RNA-expressed 
control because of its stability and constant expression in 
Salmo salar organs, in both the presence and the absence 
of the ISA virus [11,14-16]. Additionally, because of its 
relevance, ELFla has also been recommended as a house- 
keeping gene for gene expression analyses of other 
salmonid fish-related pathogens, such as Piscirickettsia 
salmonis [17], SAV [18], and VHSV [4]. 

Although an ELFla-driven assay has been successfully 
used for ISAV detection in Salmo salar samples [11,14], it 
is also known to have a highly variable performance when 
used for the diagnosis of other reared salmonid species, 
particularly Oncorhynchus kisutch and Oncorhynchus 
mykiss, which together with Salmo salar, constitute the 
economically important species cultivated in Chile. For 
this reason, we considered it imperative to analyze the 
sequence variability associated with a presumptive poly- 
morphism in the ELFla coding sequence of the three sal- 
monid species cultivated in Chile; we would then be able 
to design a transversal RT-qPCR assay that would allow us 
to use the same parameters to accurately determine the 
quantity and quality of the samples, regardless of the sal- 
monid species tested and thereby reliably diagnose RNA 
pathogens. 



Results and discussion 

Sequence variability evaluation 

Because salmon production in Chile relies on two add- 
itional species that are different from Salmo salar, we de- 
cided to compare the potential nucleotide variability in the 
RT-qPCR amplifiable target region of ELFla [11,14] in all 
three salmonid fish species reared in Chile, as variability in 
this region could be causing the species-dependent per- 
formance of the assay. cDNAs from 20 tissue samples 
from Oncorhynchus mykiss and Oncorhynchus kisutch 
were amplified, and the PCR products were sequenced 
directly to avoid any sequence selection bias that might 
result from a plasmid cloning process. The electrophe- 
rogram profiles (data not shown) consistently showed 
double peaks at four positions on the amplified region: 
one at the forward primer annealing position, one at the 
probe annealing position, and two at the reverse primer 
annealing position (Figure 1). It is particularly relevant 
that one of the two mismatches detected by the reverse 
primer was situated at the 3' end, which could seriously 
interfere with the efficiency of the amplification process 
[19]. According to previous studies, the presence of double 
peaks in the electropherograms could be explained by the 
presence of heterozygosity, SNPs or sequence insertions 
[20-22] . Although not the main focus of this work, the in- 
formation obtained here allowed us to design novel trans- 
versal primers and probes that are able to specifically 
amplify and detect selected ELFla coding sequences in 
the three fish species reared in Chile with the same degree 
of accuracy. 

To determine whether differences detected in the 
amplification could be overcome to generate a transver- 
sal RT-qPCR reaction, multiple alignments were made 
with the sequences derived from this study, as well as 
with those available in the GenBank database. Figure 1 
shows that using the BioEdit software [24], it was pos- 
sible to identify two types of nucleotide variations where 
we could detect double peaks in the electropherograms. 
First, an "intragenus conserved variation" consisted of a 
cytidine residue in the forward primer annealing pos- 
ition for all species of the Salmo genus analyzed, which 
was changed to a Thymidine residue for the two species 
of the Oncorhynchus genus analyzed (blue arrow in 
Figure 1). Additionally, we identified an "intragenus non- 
conserved variation," which did not correlate with specific 
changes in each genus (red arrow in Figure 1). 

Design of a new transversal assay 

Considering the nucleotide variations previously shown, 
we were in a position to design a new assay to detect 
ELFla in all salmonid fish species cultivated in Chile. It is 
well-known that nucleotide variations in the annealing po- 
sitions of primers and probes affect the sensitivity and effi- 
ciency of this type of reaction, as they can lead to putative 
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Forward Primer 

cm 



Salmo solar 
Salmo trutta 
Salmo marmoratus 
Oncorhynchus tshawytscha 
Oneorhynchus mykiss 
Oncorhynchus mykiss (1) 
Oncorhynchus mykiss (2) 
Oncorhynchus mykiss (3) 
Oncorhynchus kisutch (4} 
Oncorhynchus kisuteh (5) 
Oncorhynchus kisutch (6) 



Probe Reverse Primer 




E 6 i c c g 
GGGCCG 
GGGCCG 
GGGCCG 
GGGCCG 
GGGCCG 
GGGCCG 
GGGCCG 
GGGCCG 
GGGCCG 
GGGCCG 
GGGCCG 




Figure 1 Analysis of nucleotide variation between different salmonid species. A multiple alignment was made using ClustalW [23] using 
the Bioedit software [24], Top, primers and probe used in the current assay (ELFIa) [14]. The first 5 species sequences are available in the 
GenBank database; S. salar [BT060430], 5. marmoratus [FU853442], S. trutta [FF406271], 0. tshawytscha [FJ890356] ar^d 0. mykiss [AF498320]. 
Samples labeled with numbers from 1 to 6 were representative of sequences obtained from Oncorhynchus samples in this work. Blue arrow 
shows the Intragenus conserved variation and red arrow show the Intragenus non-conserved variations. 



underestimations of the actual amount of template 
present in a sample. Although it may seem trivial, mis- 
matches between primers and probes on a given tem- 
plate actually constitutes an important issue to consider 
in the generation of a highly specific and reproducible 
detection assay [19,25]. We chose to use degenerate nu- 
cleotides in the probes and primers involving internal 
nucleotide variations as a testing strategy. Additionally, 
we moved the reverse primer one nucleotide down- 
stream from the probe, taking into account that a 
variation in the 3 ' end could dramatically affect the per- 
formance of the assay. Table 1 shows the three assay 
conditions evaluated and compared in this work. Assay 
N°l (ELFla) is the assay currently in use [14]. Assay N°2 
(ELFla GIM-1), uses the same sequences as assay N°l, 
but incorporates degeneracy in the variable positions. 
Assay N°3 (ELFla GIJVI-2) employs different primers to 
eliminate the 3' end degeneration in the original reverse 
primer. The same probe was used for assays N°2 and N°3, 
which differ from each other only in terms of the primers 
used. 



Evaluation of the alternative new assays 

Initially, performances of the two novel alternative assays 
(N°2 and N°3) were compared against the classical one 
(N°l) using templates obtained from S. salar, O. kisutch 
and O. mykiss tissue samples. For this comparison, the 
detecting probe was used without an MGB or LNA 
system [13]. Under these conditions, the probe displayed 
a Tm value of 49.9°C, which was calculated using 
OligoCalc [26] and which suggested low stability in the 
annealing process of the probe and template. This 
allowed us to evaluate the performance of the assay 
based on the condition that any mismatch could affect 
the results considerably, thus making it easier to detect 
any improvement. 

The results of the three assays in Table 2 clearly show 
that assay N°3 performed better than assays N°l and N°2 
and was therefore selected for further evaluation. Assay 
N°3 had lower Ct values and higher relative fluorescence 
than the other two assays for all samples evaluated. As 
expected, the Ct value differences between assays N°l 
and N°3 were higher between samples from O. mykiss 



Table 1 Alternative assays used to compare efficiency and transversality of detection 



Assay N° 


Assay name 


Name 


Sequence (5' - 3') 


Reference 


1 




FLFIa For 


CCCGTGCAGGACGTTTACAAA 


[14] 


1 


FLFlo 


FLFIa Rev 


CACACGGGCCACAGGTAGA 


[14] 


1 




FAM-FLFla 


FAM- ATGGGCGGTATTGGAAG 


[14] 


2 




FLFIa GIM-l-For 


CGCCTCGAGGAYGTYTACAAA 


This work 


2 


FLFIa GIM-1 


ELFla GIM-1 Rev 


CACACGGCCCACRGGTAGW 


This work 


2 




FAM- ELFla GIM 


FAM-ATCGGYGGTAT 4-T-fG-fG + A + AC-BHQ 


This work 


3 




ELFla GIM-2 For 


GCGCCTCCAGGAYGmACAA 


This work 


3 


FLFIa GIM-2 


ELF a GIM-2 Rev 


GCAGACGGCCGAGRGGTAC 


This work 


3 




FAM-FLFla GIM 


FAM-ATCGGYGGTAT 4-T4-G4-G + A + AC-BHQ 


This work 



Assay N°1 (ELFla) is the current assay. Assays N°2 (ELFla GIM-1) and N°3 (ELFla GIM-2), use alternative primers and probes designed to test the objectives of this 
study. The FAM-ELF1 a probe has MGB, and the FAM-ELF1 a-GIM probe has LNA residues, denoted by a -I- symbol before the corresponding nucleotide. 
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Table 2 Comparative performance of the assays in different species of salmonid tested 

ELF1 a ELF1aGIM-1 ELFIa GIM-2 

Ct value ARn Ct value ARn Ct value ARn 

O.mykiss 32.31 ±0.88 0.201 ±0.039 26.49 ±0.25 0.811 ±0.019 21.43 ±0.080 1.498 ±0.01 7 

O.kisutcb 31.10±0.36 0.238±0.021 27.86±0.17 0.686±0.043 22.71 ±0.109 1.353±0.019 

S.salar 28.55±0.43 0.307±0.021 2918±0,14 0.519±0.008 23.56±0.017 1.232±0.001 

The comparison was made based on the Ct values and relative fluorescence intensities (ARn). The probes used in this comparison were depleted of either MGB or 
LNA. The data are shown as a Ct value ± standard deviation and ARn ± standard deviation (n = 3). 



and O. kisutch species. This is because assay N°3 takes 
into account the nucleotide variability in the amplifica- 
tion region for this genus, thereby avoiding the mis- 
matches observed in assay N°l. In the same way, 
samples from S. salar showed slighter differences in Ct 
values between assays. 

Although both alternative assays (N°2 and N°3) con- 
sidered the existing nucleotide variability of the species 
tested, the performance of assay N°3 was notoriously 
better than that of assay N°2. We hypothesize that this 
effect may be a result of a structural interference be- 
tween the primer and the probe and not primarily of the 
template mismatches. 

The next step in the evaluation was to consider the de- 
sign of a more stable assay by increasing the Tm of the 
probes and thus increasing the annealing stability between 
the probe and the template. To achieve this objective, 
MGB was used on the probe of the ELFla assay, as 
reported previously [13] and Locked Nucleic Acid (LNA) 
was used on the probe of the ELFla GIM-2 assay. The 
comparison was made between assay N°3 (selected new 
assay) and the original assay N°l (reference assay) consid- 
ering samples derived from both tissues and cell cultures 
from different species of salmonids. Table 3 summarizes 
our results, showing that, as expected, cell culture and tis- 
sues samples from S. salar had lower ACt values between 
assays (Ct value ELFla - Ct value ELFla GIM-2), while 
cell culture and tissues samples from the Oncorhynchus 
genus had higher differences in Ct values between assays 
and therefore a higher ACt. The largest differences. 



however, were observed between samples derived from O. 
tshawytscha, where a maximal ACt value of 11.101 was 
obtained from tissue-derived samples. 

An additional analysis can be performed when we 
compare the results from Table 2 (probes without the 
enhancer Tm system) with the results from Table 3 
(probes with the enhancer Tm system). In fact, differ- 
ences between assays N°l and N°3, shown in Table 2, 
gave a ACt value of 10.88 for the O. mykiss tissue sam- 
ple, whereas in Table 3, the ACt value was diminished to 
5.27. This is a powerful demonstration that the en- 
hancers, in addition to increasing Tm values, stabilize 
the probe-template annealing and, as a result, render the 
assay more robust. 

Although only LNA was used in the probe of the 
ELFla GIM-2 assay, we would expect a similar result 
when using MGB as an enhancer. 

The evaluation of the efficiency of assay N°3 for all 
samples tested is shown in Table 3. All values for the 
RT-qPCR reactions fell within the suitable experimental 
range (90%-110%) [27]. 

Evaluation of field samples 

The final evaluation compared assay N°3 (ELFla GIM-2) 
with the original assay N°l (ELFla), using the same dis- 
tribution of LNA as a system to increase the probe Tm 
value in both assays. Field samples from S. Salar (n = 98), 
O. mykiss (n = 97) and O. kisutch (n = 82) were obtained 
from different farms located in the south of Chile, and 
parallel assays were run under the same RT-qPCR 



Table 3 Comparison of differences in Ct values in tissue and cell culture samples from different salmonid species 



Samples 


Type 


Specie 


ELFla 


ELFla GIM-2 


A Ct 


Eficiency 








Ct Value 


Ct Value 


value 


(%) 


ASK-1 


Cell culture 


S. salar 


28.1 34 ±0.1 49 


25.381 ±0.114 


2.753 


96.65 


SHK-1 


Cell culture 


S. salar 


27.1 98 ±0.204 


24.237 ± 0.078 


2.961 


92.78 


RTS11 


Cell culture 


0. mykiss 


30384 ± 0338 


25.028 ±0.1 91 


5.356 


96.52 


CHSE-214 


Cell culture 


0. tshawytscha 


32.768 ±0.421 


24.354 ± 0.262 


8414 


93.64 


GIMCP 


Tissue 


S. salar 


29.235 ±0.1 07 


25.297 ± 0.205 


3.938 


96.90 


GIM017 


Tissue 


0. mykiss 


33.467 ± 0.63 


26.1 94 ±0.1 99 


5.273 


96.16 


GIM025 


Tissue 


0. kisutch 


29.447 ±0.1 77 


24.693 ±0.1 46 


4.754 


101.24 


GIM496 


Tissue 


0. tshawytscha 


36.494 ±0.2 17 


25.393 ±0.161 


11.101 


97.22 



The ACt value was calculated (Ct value ELFla - Ct value ELFla GIIVl-2). 
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conditions (reagent mix, software setup, and data analysis). 
Four (out of 82) samples from O. kisutch were withdrawn 
from the analysis because they could interfere with the 
correct analysis because they displayed ACt values over 9, 
and one of them had a ACt value of 16. This result sug- 
gested that this might be due to several mismatches be- 
tween the template and the primers and the probe of the 
ELFla assay. 

An integral analysis of the data shows that, similar to 
previous results, S. solar samples displayed smaller differ- 
ences between assays N°3 and N°l, which was expected, as 
the latter was designed to specifically fit S. salar se- 
quences. Samples from O. mykiss and O. kisutch displayed 
higher differences in Ct values, as seen in Figure 2. There- 
fore, ACt values depend on the salmonid species analyzed. 
Furthermore, no significant differences in ACt values were 
observed between the O. mykiss and O. kisutch samples 
analyzed. 

Figure 3 shows the Ct values displayed in each assay 
from field samples of the three salmonid species. As 
expected, the new assay N°3 (ELFla GIM-2) displayed 
no significant differences among the Ct values obtained 
(Figure 3B). When the original assay N°l (ELFla) was 
used, significant differences were observed between Ct 
values of samples from the S. salar and Oncorhynchus 
species (Figure 3A). 

These differences in Ct values reveal the importance 
in diagnoses of an adequate endogenous control. Indeed, 
using assay N°l, samples derived from S. salar display Ct 
values within a valid range for diagnostic purposes, 
whereas samples from the Oncorhynchus group display 
higher Ct values, some of which were beyond the valid 





6,0-1 




5,5- 




5,0- 




4,5- 




4,0- 




3,5- 


OJ 




D 






3,0- 


> 






2,5- 


u 




2,0- 




1,5- 




1,0- 




0,5- 




0,0- 



S. salar O. mykiss O. kisutch 

Figure 2 ACt values from the ELFla GIM-2 and ELFla assays In 
field samples. The ACt value was calculated (Ct value ELFla - Ct 
value ELFla GIM-2). In the plot, error bars are standard deviations 
from the average of all data. S. satar is significantly different from the 
0. kisutch and 0. mykiss species, according to the one-way ANOVA 
following by a Tukey's test (p < 0.05). 



25 



20- 



> 

4-' 

U 



S.salar 



O.mykiss 



O.kisutct) 



B 



25-, 



20- 



0) 

0! 
> 



15 



10 



S.salar 



O.mykiss 



O.kisutch 



Figure 3 Distribution of Ct values of field samples. ELF1 a assay 
(A) and the ELFla GIM-2 assay (B). The ELFla GIM-2 assay does not 
show significant differences between the three salmonid species, 
whereas the ELFla assay shows a significant difference between the 
S. salar and Oncorhynchus species samples (x). Statistical analysis was 
performed using one-way ANOVA and Tukey's test (p < 0.05). 



range for diagnoses. Thus, incorrect diagnoses may 
occur as a result of a weak detection assay design; how- 
ever, this effect is corrected in the new assay. 

Conclusions 

Here we identified nucleotide variability in the amplified 
region of ELFla, which causes highly variable results when 
the ELFla assay is used as an endogenous control for 
Oncorhynchus species. Using this information, we have 
designed and improved a single reliable and efficient assay 
(ELFla GIM-2), driven to an endogenous cellular quality 
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control, to be used in RT-qPCR diagnostics of pathogenic 
agents in tissue samples from the salmonid species ana- 
lyzed in this work. 

Methods 

Samples 

Tissue samples analyzed were obtained from the kidneys 
of Salmo salar, Oncorhynchus mykiss, Oncorhynchus 
kisutch and Oncorhynchus tshawytscha according to offi- 
cial procedures established by Sernapesca, the Chilean 
government institution in charge of animal health sur- 
veillance in aquaculture and fisheries. Salmonid cell lines 
ASK-1, Atlantic Salmon Kidney (ATCC; CRL-2747) [28], 
SHK-1, Atlantic Salmon Head Kidney [29], CHSE-214, 
Chinook Salmon Embryo (ATCC; CRL-1681) [30] and 
RTSll, rainbow trout monocytes/macrophages cell line 
[31], were used as verification controls. 

Ethical statement 

Tissue samples from fish were obtained from the sur- 
veillance program for fish disease in Chile. Fish were not 
killed for the purpose of this study. All sampling was 
performed according the regulations of Sernapesca 
(Chilean government institution in charge of fish health) 
and carried out in strict compliance with the recommen- 
dations of Chapter 7.4 of the Aquatic Animal Health 
Code of the World Organization for Animal Health. 
Every effort was made to minimize animal suffering in 
all procedures. 

Total RNA Isolation 

Tissue samples were preserved in RNAlater (Ambion) to 
be transported from the fish farms to the laboratory. 
Total RNA from the tissue samples and cell culture sam- 
ples were isolated using the RNeasy Mini Kit (Qiagen), 
according to the standard protocols recommended by 
the supplier. Prior to the isolation of RNA, tissue sam- 
ples were homogenized using a MagnaLyser device 
(Roche) for 30 sec at 6500 rpm. RNA samples were 
stored at -80°C until use. 

cDNA synthesis 

RNA samples for sequencing were reverse transcribed 
using the Superscript III kit (Invitrogen), according to 
the protocols recommended by the supplier. Briefly, ap- 
proximately 1 |ig of total RNA was added to the mix, 
which contained 0.5 ^iL of dNTPs (10 mM), 2 |iL of 5 X 
First Strand Buffer, 0.5 ^iL of DTT (0.1 M), 0.5 |iL of 
Superscript III reverse transcriptase enzyme, and 1 i^L 
of Random Primers (100 |iM) (Fermentas) in a final vol- 
ume of 10 i^L. Reactions were incubated at 25°C for 
10 min, 50°C for 1 h, and 75°C for 15 min. The final 
cDNA was diluted 3-fold before being used as a tem- 
plate in PGR. 



PCR and sequencing analyses 

All PCR reactions for sequencing purposes were per- 
formed using external primers ELFla Ext For 5'-ATG 
GGC TGG TTC AAG GGA TG-3' and ELFla Ext Rev 
5 -GGT GGT GCA TCT CCA GAG AC-3'. PCR was 
performed using the Go Taq Flexi Polymerase Kit 
(Promega). Each reaction had a final volume of 25 i^L and 
consisted of 5 |iL of 5X Colorless Flexi Buffer, 0.5 i^L of 
dNTP (10 mM), 3 ^iL of MgCLj (25 mM), 0.2 |iL of Go 
Taq DNA Polymerase, each primer at a final concentration 
of 400 nM, 11.3 |iL of nuclease-free water (Invitrogen), 
and 4 i^L of the diluted cDNA. The PCR program used 
consisted of an initial denaturation for 2 min at 95°C, 
followed by 40 cycles of denaturation for 30 sec at 94°C, 
primer annealing for 30 sec at 60°C, and elongation for 
1 min at 72°C. The final elongation step was run for 5 min 
at 72°C. The resulting amplicons were resolved by agarose 
gel electrophoresis, and the target bands were purified 
from the gel using the E.Z.N.A Gel Extraction Kit (Omega 
Biotek), according to the protocols recommended by the 
supplier. The purified PCR products were quantified and 
then sequenced at Macrogen Inc. (Korea) with the same 
primers used in the PCR reaction. Sequences were inter- 
preted using the Bioedit Software [24] and the multiple 
alignments were performed via ClustalW [23]. 

Real-time RT-qPCR 

Real-time PCR amplifications were performed using the 
Super Script III Platinum One-Step Quantitative RT-PCR 
System Kit (Invitrogen) and the StepOnePlus Real-time 
PCR System thermocycler (Applied Biosystems). Primer 
sets and assay names are listed in Table 1. Each reaction 
was carried out in a final volume of 20 \iL. Primers were 
at a final concentration of 1 [^M, and hydrolysis probes 
were at a final concentration of 0.3 [iM. The reactions 
contained the passive reference dye, ROX. One cycle of re- 
verse transcription was run for 15 min at 50°C, 1 cycle of 
denaturation for 2 min at 95°C, and 45 cycles of denatur- 
ation for 15 sec at 95°C, followed by annealing and elong- 
ation for 1 min at 60°C. Real-time PCR efficiencies were 
calculated from the slope according to the established 
equations = 10 (-i/'=i°p<=) [32]. 

Field sample evaluation 

Samples consisted of organ pools (heart, gill and kidney) 
from 98 specimens of S. salar, 97 of O. kisutch and 82 of 
O. mykiss. Samples were fixed with RNAlater in a tissue: 
fbcator ratio of 1:10 and stored at 6-8°C prior to analysis. 
The fixed samples were homogenized using Precellys de- 
vice (Bertin Technologies). The homogenized samples 
were centrifuged by 2 min. at 13.000 r.p.m (pico 17, 
Thermo). The supernatants were purified for RNA using 
a high purity viral nucleic acid kit (Roche) according to 
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the manufacturer's instructions. The purified RNA was an- 
alyzed using the RT-qPCR protocols described previously. 

Statistical analysis 

The statistical analyses were performed using OriginPro 
V8.5 software. Comparisons between Ct values or ACt 
values were performed using one way ANOVA followed by 
a Tukey's tests (P < 0.05), to denote significant differences. 

Abbreviations 

ELFla: Elongation factor 1 alpha; PCR: Polymerase chain reaction; 
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